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Abstract 

We investigate ground-state and excitation spectrum of a system of 
non-relativistic bosons in one-dimension interacting through repulsive, 
two-body contact interactions in a self-consistent Gaussian mean-field 
approximation. The method consists in writing the variationally de- 
termined density operator as the most general Gaussian functional of 
the quantized field operators. There are mainly two advantages in 
working with one-dimension. First, the existence of an exact solution 
for the ground-state and excitation energies. Second, neither in the 
perturbative results nor in the Gaussian approximation itself we do 
not have to deal with the three-dimensional patologies of the contact 
interaction . So that this scheme provides a clear comparison between 
these three different results. 



PACS numbers : 02.30.Mv, 05.30.-d, 05.30.Jp, 67.40.Db 

1 Introduction 



The non-relativistic interacting Bose gas is certainly among the most inter- 
esting problems in many-body physics. This system has been extensively 
investigated since fifties motived strongly by the superfluid properties of he- 
lium [|l|. The recent successful achiement of Bose- Einstein condensation of 
alkalis vapors has renewed a great interest in theoretical studies of Bose 
gases Ij. 

For dilute gases the macroscopic condensate function is given by the 
Gross- Pit aevskii equation. The effects of noncondensate atoms can be cal- 
culated in the loop expansion in power of y/pa^, where p is the density and 
a is the S'-wave scattering [^. This perturbative method tends to be less 
reliable in the circumstance when higher trap densities are involved. On the 
other hand, non-perturbative approximation schemes have been developed 
in connection with the problem of self-interacting, relativistic Bose fields, 
which became relevant for inflationary models |^, and in the analysis of the 
transient phenomena in the collision of the complex nuclear system [^ . This 
problem has also been used as a testing ground for non-perturbative meth- 
ods which could then be applied to more complicated systems of interacting 



fields. In this context, the Gaussian variational approach [0 has received 
considerable attention, also in view of its relation to extended mean field 
methods traditionally employed in many-body physics g, |^. 

The main objective of this work is to compare the standard perturbation 
results with the unperturbative ones obtained with the Gaussian approxi- 
mation in the case of a system of non-relativistic bosons in one dimension 
interacting through repulsive, two-body contact interactions. The use of 
an one- dimensional problem gives the advantage of dealing with a repulsive 
contact interaction without getting into problem due to divergencies ll^, |ll|]. 
Also the existence of an exact solution in this case show us clearly when the 
perturbation theory breaks down and how the upperbound non-perturbative 
results stands in this region. 

The Gaussian mean-field approximation which we will employ here is 
based in terms of a time-dependent projection approach developed earher 
for the nonrelativistic nuclear many-body dynamics by Nemes and de Toledo 
Piza [Q. This scheme allows one to formulate a mean-field expansion for 
the dynamics of the two-point correlation function from which one recov- 
ers the results of the Gaussian mean-field approximations in lowest order, 
i.e., this permits to include and to evaluate higher dynamical corrections 
effects to the simplest Gaussian mean-field approximation. Moreover, the 
expansion is energy- conserving (for closed system) to all orders [0]. The 



method was recently applied in the context of the field-theoretical models, 
namely, to treat an uniform relativistic (H-l)-dimensional self-interacting 
boson system described by the A0^ theory |p, to treat an uniform relativis- 
tic (H-l)-dimensional self-interacting fermion system described by the chiral 



Gross- Neveu model [Q, and to treat an uniform (3+l)-dimensional relativis- 



tic interacting boson-fermion system described by the Plasma Scalar model 



Although we developed a temperature dependent formalism in this paper 
we will concentrate on the zero temperature case and calculate in Gaussian 
mean-field approximation the physical quantities related to the optimized 
free energy such as ground state energy and sound velocity to compare with 
pertubative and exact results given in literature [|^, |T^. In Section 2, we 
define the problem and show its exact solution. Section 3 reviews the approx- 
imation used here and discusses the equihbrium solutions for the resulting 
variational equations. Through a truncated version of the Gaussian varia- 
tional equation we obtain the perturbative results in Section 4. In Section 5 



we show the numerical results obtaneid from the Gaussian and perturbative 
approximations and compare these with the exact calculations. Section 6 
discusses the range of validity of these theoretical methods and main conclu- 
sions. 



2 The problem and the exact solution 

We consider an uniform and isotropic system of non-relativistic, interacting, 
spinless bosons in one-dimension described by the Hamiltonian, 

^=-E(^)+2cE5(x.-x,) (1) 

with periodic boundary condition in region < Xi < L, where 2c > is the 
amplitude of the repulsive delta function potential. The only dimensionless 
intensive variable in the theory is 

c 

7 = - , 
P 

where p = N/L is the density system. 

Lieb and Liniger in Ref. [jl6[ have shown that the exact ground state 



energy can be obtained by solving an inhomogeneous Fredholm equation of 
the second kind. So that in the thermodynamic limit (A^, L ^ oo such that 
p = fixed constant) we can write the ground state energy per particle as 

^° - ^ ^"^ f{k)edk , (2) 



A^ pJ-K 

where f{k) is the solution of 



with the condition 

''^f{k)dk = p . (4) 

For computational purposes we can change variables as follows 



k = Kx ; c = K\ ; f{Kx) = g{x) 



in terms of which the Eqs.(2), (3) and (4) become, respectively, 

-J^ = P'^J^ J_^9{x)x^dx (5) 

l + 2A/;^j^i!^rf. = 2.9fe) (6) 

7 / g{x)dx = A . (7) 



-1 

For a fixed A we can calculate g{x) having A(7) and with it follows ^{'^) and 
-^(7) EH- Once we have the ground state energy it is possible to calculate 
the chemical potential 

From /i we can calculate the sound velocity derived from the macroscopic 
compressibility 



Lieb, in another work [|17l, also has shown that this theory exhibit two types 



of excitations. The energy of the first type of excitation is given by 



e-^{p) = -^ + q^ + 2 kJ{k)dk , (10) 

J —K 



where /i is the chemical potential and 



p = q+ f J{k)dk , (11) 



-K 

while J{k) is again a solution of a Fredholm equation 



271 J{k) =2c f - — ^T— — TT - TT + 2 tan-i(g - k) . (12) 

J-K & + (k — rV 



-K c^ + [k — r)^ 



The energy of the second type of excitation have the same structure and is 
such that 

e2(p) =/x-g2 + 2 / kG{k)dk (13) 

p=-q+ f G{k)dk (14) 



2nG{k) = 2c r , ^y^^"' ,, + TT - 2tan-^(g - k) . (15) 

J-K C' + [k — ry 

To solve these equations we use the same procedure used in the ground state 
energy. Here q plays the same role as A in the other case. 



3 Gaussian approximation 

In second quantization, rewriting the Hamiltonian (1), in momentum repre- 
sentation with periodic boundary conditions in region L, follows 

H = ^ e{k)a\'ar:+ — \^ a\-^ jx]-. jxcac , (16) 

k kik2q 

where e{k) = k'^ is the free particle kinetic energy, and the contact repulsive 
(c > 0) interaction potential between a pair of particles is 2c 6{x — x'). The 
at and a^ are boson creation and annihilation operators, which satisfy the 
standard boson commutation relations at equal times. 

3.1 Projection technique and approximation scheme 



In this section we introduce the time-dependent projection technique |]T2 
which will permit to obtain closed approximations to the dynamics of a sys- 
tem. It has been developed earlier in the context of nonrelativistic nuclear 
many-body dynamics and was recently applied in the quantum-field theoreti- 
cal context [^, |T^, . It allows for the formulation of a mean-field expansion 
for the dynamics of the two-point correlation function from which one recov- 
ers the results of the Gaussian mean-field approximations in lowest order. If 



carried to higher orders it allows for the inclusion and evaluation of higher 
dynamical correlation corrections to the simplest mean-field approximation. 
We begin by describing the quantum state of the system in Heisenberg 
picture in terms of a many-body density operator JF, a time independent, 
non-negative, Hermitian operator with unit trace. T may in general involve 
statistical mixtures, in particular in the boson number operator, in the spirit 
of a grand canonical description. Our implementation of the Gaussian ap- 
proximation will consist in decomposing the full density T as 

T = T^{f)^T\t) , (17) 

where T{){i) is a Gaussian ansatz which achieves a Hartree-Fock factorization 
of traces involving more than two field operators. The most general hermitian 
Gaussian density T^if) have the form of a exponential of linear and bilinear 
expressions in the fields normalized to unit trace []TB|, 0. In the momentum 
basis, it reads 



T. 



exp (E(fc;,fe) {^^^.sA.'^k, + [^k.^AA^, + h-c-) + (^fcl4i + ^-C' 



Tr {exp (E(,-;,,-;) {A^,sA\''^. + [^k.sAA, + ^•^•) + fe4i + h-C-)])} 

(18) 



% = ^Ih + yV^_k 






where the matrix A^ ^ is hermitian. The parameters in Eq.(18) are fixed by 
requiring that mean values in Tq of expressions that are bilinear in the fields 
reproduce the corresponding T averages [see Eqs.(21) and (22) below]. 

The uniformity and isotropy assumptions we make allow us to restrict 
this general ansatz so that A^_^ ^- is diagonal, B^^ ^- vanishes unless k\ = — ^2 
and both of these matrices and the C^- depend only on the magnitudes of 
the momentum vectors. Finally, by a general canonical transformation of 
the Bogolyubov type, the time- dependent matrix JFq can be diagonalized [§]. 
In this case, the diagonalization is achieved by defining transformed boson 
operators as 



where 



h = «fc - r* 



Tfc = (ag) . 

We also have used the isotropy of the uniform system to make the c-number 
transformation parameters Xk, Hk and Tk dependent only on the magnitude 
of k. In order for this transformation to be canonical we have still to impose 
on the Xk and yk the usual normalization condition 

\xk\^ - \yk\^ = I . (19) 

In terms of the Bogolyubov quasi-boson operators, the trace-normalized 
Gaussian density operator is now written explicitly as 

^0 = n -, -r^ ■ 20 

Straightforward calculation shows that 

Trivl^Vk-:.J'o) = i^kAki - h) (21) 

so that the Uk are positive quantities corresponding to mean occupation num- 
bers of the r^-bosons. One also finds that 

Tr[{xla^, + y^ot )"^o] = Tr[{Vk + ^fc)"-^o] = {AkT (22) 

with Ak = xlTk + Vk^k so that non vanishing values of the Tk correspond 
to coherent condensates of unshifted Bogolyubov transformed bosons. We 
again invoke the uniformity of the system to impose 

Tfc = SkflTo (23) 

in the calculations to follow. 

The "remainder" density J-''{t), defined by Eq.(17), is a traceless, pure 
correlation density. As already remarked, a crucial point to observe is that 
jFo(t) can be written as a time-dependent projection of JF, i.e.. 



J^^{t) = V{t)T with V{t)V{t) = V{t) . 



(24) 



In order to completely define this projector we require further that it satisfies 
iTo{t) = [V{t), L\T= [^o(t), H] + V(f) [H, T\ , (25) 



where C is the Liouvillian defined as 



C- = [H, 



(26) 



H being the Hamiltonian of the field. Eq.(25) is just the Heisenberg picture 
counterpart of the condition V{t)J-' = which has been used to define V{t) 
in the Schrodinger picture |jl3|]. It is possible to show that conditions (24) 
and (25) make V{t) unique and to obtain an explicit form for this object in 
terms of the quasi-boson operators and of the natural orbital occupations see 
Refs. i, in, ill. 



The existence of the projector V{t) allows one to obtain an equation 
relating the correlation part J^'it) to the Gaussian part J-'o(t) of the full 
density. This can be immediately obtained from Eqs.(17), (24) and (25) and 
reads 



(tdt + V{t)C) J'' it) = (J - V{t)) CJ^oit) 



(27) 



This equation has the formal solution 



^'(t) = g(t, o)r (0) - 1 f dt' g{t, f) (J - p(t')) cj^oit') 

Jo 



(28) 



where the first term accounts for initial correlations possibly contained in J-". 
The object G(t,t') is the time-ordered Green's function 



Q{t,t') = T f exp i j 



I I dTV{T)C 

t' 



(29) 



We see thus that J^'{t), and therefore also JF [see Eq.(17)], can be formally 
expressed in terms of J^oit') (for t' < t) and of initial correlations JF'(O). 
This allows us to express also the dynamics of the system as functionals of 
J^o{t') and of the initial correlations. Since, on the other hand, the reduced 
density J-'o(t') is expressed in terms of the one-boson densities alone, we see 
that the resulting equations are essentially closed equations. Note, however, 
that the complicated time dependence of the field operators is explicitly 
probed through the memory effects present in the expression (28) for J^'{t). 
Approximations are therefore needed for the actual evaluation of this object. 
A systematic expansion scheme for the memory effects has been discussed 



in Refs. p, |1^, [1^, |T5[. The lowest order correlation corrections to the 
pure mean field approximation, in which JF' is simply ignored, correspond 
to replacing the full Heisenberg time-evolution of operators occurring in the 
collision integrals by a mean-field evolution governed by 

H, = V\t)H . 

Consistently with this approximation, C is replaced in (25) and (26) by £o " = 
[ifo, ■ ]. In this way correlation effects are treated to second order in H in 
the resulting collision integrals. 

An important feature of this scheme (which holds also for higher orders 
of the expansion |[1^) is that the mean energy is conserved, namely 

where 

{H) = Tr HTo{t) + Tr HT'{t) . 

In the following sections we apply the general expressions obtained in 
above to treat a uniform boson system described by the Hamiltonian given 
in Eq.(16). We will consider only the lowest (mean-field) approximation, 
corresponding to J-''{t) = 0. Collisional correlations will be treated elsewhere. 

3.2 Temperature time-dependent Gaussian treatment 

First, it is important to note that the truncated density JFq in general breaks 
the global gauge symmetry of H, which is responsible for the conservation of 
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the number of a-bosons. It is possible to verify that assumption calculating 
the dispersion 



k,k' k 

Calculating the traces above we have 

(Ar2)_(iV)2 = 2\ro\'[\xo\'uo + \yo\\l + uo)]-2rfxoy*o{l + 2uo)- 

-2rlyoxl{l + 2z/o) + iFoP + Y.[\^k\^^k + (1 + UkMl^] + 

k 



Furthermore, mean values of many-boson operators taken with respect to JFq 
will contain no irreducible many-body parts, so that the replacement of JF 
by J-'o amounts to a mean field approximation. The states described by J-'o 
have therefore to be interpreted as "intrinsic" mean field states. 

In order to develop a finite temperature treatment within this framework, 
we look for variational extrema of an approximation Q to the grand potential 
written in terms of the Gaussian density operator jFg as 

n = Tr[J^o{H-i2N + KT\nJ^o)] (30) 

with 

it 

where K is the Boltzmann constant and the Lagrange multipliers /i and T 
will play the role of chemical potential and temperature respectively. The 
replacement of the full density operator JF by jFg implies not only truncation 
of correlation energies but also replacement of the full entropy by the mean 
field entropy 5*0 = —KTr[J^o\nJ^o]. Variations will be taken with respect to 
the parameters which determine the Gaussian density J-'o, and the value of /i 
is fixed by choosing the particle density of the system. A detailed justification 
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of this general procedure has been given [|T0|, |T9], who have shown that it 
corresponds to the optimal determination (in the variational sense) of the 
grand partition function (hence of the grand potential) for the problem when 
one restricts oneself to trial density operators of the form given in Eq.(20). 
We have to emphasize that the use of this functional Q does not assure that 
we will have the best variational results for other observables. See discussion 
in Sees. 5 and 6. 

An important simplification occurs in the case of stationary problems. In 
this case Xk, yu and Fq can be taken to be real and we can use a simple para- 
metric representation that automatically satisfies the canonicity condition. 
It reads 



Xk = cosh (jfc , yk = sinh cxfc . 



(31) 



It is then straightforward to evaluate the traces involved in Q to obtain 



n 



Y. e(^) -^^ 



Acvr 



[1 + 2z/,fc) cosh2cTjfc — 1 



•^X:(l + 2^fc)sinh2afc + |^^ 



^r^ + Y 



(1 + 2z/fc) cosh2crA,. — 1 



+ ^<T.i^ + '^^k)smh2aA (32) 

-fCT^[(l + z/fc)ln(l + z/fc) -z/felnz/fc] . 

k 

In a similar way the number constraint Tr[NJ-'Q] = (N) evaluates to 

(iv) = rg + x:[^'^''^^7^'"'"'1 . (33) 

k L 



3.3 Equilibrium solutions 

Equations determining the form of the Gaussian density J-'o appropriate for 
thermal equilibrium are in general derived by requiring Q, given by (32), to 
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be stationary under arbitrary variations of Fq, (Jk and z/^. Variation with 
respect to Fq gives the gap equation 

Fo|yF2-2/i-4cA + 8cfi| = , (34) 

where we define in the thermodynamic hmit, J2k ~^ ^ I-^ dk, the quantities 
A, B and C given by 



1 f°° 

— (1 + 2z/fc) sinh 2(7^ rf/c 

m J-oo 



All J-OO 

1 /"^ 
B = — [(1 + 2z/fc) cosh2afe - 1] rfA; 
47r J-oo 

1 /""^ 
C = — e{k)[{l + 2uk)cofi\i2ak-l]dk . 

4:71 J-OO 

The gap equation, Eq.(34), admits a solution with a non- vanishing value 
of Fq (condensed phase) obtained by requiring that the expression in curly 
brackets vanishes. This condensed phase solution involves the number con- 
straint, Eq.(33), in addition to the values of Uk and 0"^, which are determined 
by the remaining variational conditions on Q. In order to simplify the al- 
gebraic work involved in the study of this class of solutions, it is sometimes 
convenient to use the Eq.(33) to eliminate Fq from Q, which then becomes 

n = F-ij{N) . 

This identifies the free energy F as 



F = LC-2cp{A-B)L + c{p'^ + A'^-B^)L + 2cABL 

(35) 

LKT /•°° 
— / [(1 + z/fc) ln(l + z/fc) - t'fc In z/fc] dk . 

For condensed phase, the chemical potential follows from Eqs.(33) and (34) 

13 



fi = 2c{p-A + B) . 



(36) 



Extremizing F by setting derivatives with respect to Ck and Uk equal to zero 
one gets 



and 



where 



tanh 2o"r 



2c[p -B-A] 
e{k) + 2c[p -B + A] 



T^k 



{exp[^fE/KT] - 1} 



(37) 



(38) 



eg{k) = v^ = ^Je{k)^ + Ace{k)[p - B + A] + 16c2[p - B]A (39) 

corresponds to the excitation spectrum in our approximation. This result 
can be confirmed through an RPA calculation Ref. [^]. Finally, A and B 
can be found from gap equations (37-38), for given values of p and T, by 
solving numerically the system of equations 



A = — 



1 

4tt 



( 2c[p -B -A] 



1 + 



{exp[VA/KT] - 1}_ 



dk 



B 



1 /-^ ( e{k) + 2c[p - B + A] 

AtT 7-00 1 



A 



{exp[VA/KT] - 1}_ 



(40) 
l\dk . 



Therefore, in terms of the values of A and B, we can calculate C and the 
thermodynamic functions also numerically. 

On the other hand, the gap equation (34) may also admit, besides the 
condensed phase, the trivial solution Fq = 0, which corresponds to a non- 
condensed phase. For this phase, we extremize fl, given in Eq.(32), by setting 
derivatives with respect to ak and Uk equal to zero one gets 



tanh 2a k 



-2cA 



e{k) - p + 2cB 



(41) 
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J^k = ^ , (42) 

{exp[VA/KT]-l} ^ ' 



and 



where 

eg{k) = VA = e{k) - ^i + Acp . (43) 

For non-condensed phase (Fq = 0) the only possible solution, from gap equa- 
tions (41-42), is A = following that cTk = 0. Therefore, from Eqs.(33) and 
(42-43), we have the solution 

^~27rJ-ooexp[{k^-l2 + Acp)/KT]-l ' ^ ^ 

which serves to determine fi. That is the standard Hartree-Fock result which 
corresponds to a shift in the chemical potential when compared with the 
ideal Bose gas [^ . In this phase, we show in Ref. ^^ that when T ^ we 



have n — > 4cp. 

Particularly in this paper we are interested on properties for T = 0, since 
in this case we can compare our results with perturbative and exact solutions 



existent in the literature [^, ^7j. When T = 0, the ground state energy is 



then given by 




F{T = 0) 


E 


N 


N 



--2c(A-B) + -(p^ + A^-B^) + 2-AB . (45) 
P P P 

We still have to decide what phase shall we use when T = 0. To do so, 
let us examine the chemical potencial p since from Eq.(36). When A > 
follows from Eq.(36) that yUroT^o < 4cp = yUro=o and it is easy to check that 
for a given p at T = there are always A and B that are solutions of the 
system (40). This does not happen in the case of non-condensed phase at 
T = [|1^. Therefore, we conclude that the stable phase at T = for any p 
is the one with Fq 7^ 0. 

Finally, in this calculation A,B and C are finite meaning that no renor- 
malization is needed. This is a quite different result from the three-dimensional 



case Ref. m 20 
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4 Independent 77-bosons - Perturbative results 



In Ref. [|T0| it is shown that from a convenient truncation of the Gaussian ap- 
proximation in the three-dimension problem, it is possible to obtain the usual 
perturbative results in the parameter 2c. Here we will follow the discussion 
of Ref.[|l|,|ll. 

This truncation of the Gaussian approximation consists in neglecting all 
terms representing interactions between ?7-bosons, what amounts to dropping 
all double integrals in Eq.(35). The variational conditions on ak and z/^ 
appear then as 



tanh 2o"i 



2cp 



k^ + 2cp 



(46) 



and 



T^k 



{expyk" + icpk'^/KT] - 1} 



(47) 



which gives us the usual phonon excitation spectrum 



e(p) = Jp'^ + Acpp^ 



(48) 



When we using these results for calculating F, again stressing that there is 
no need for regularization here, we obtain 



E 2 L 



V7' 



(49) 



Using Eq.(8) and Eq.(9) for the chemical potential and the sound velocity 

we get 



.^^M-^y 



Vs = 2p 



1 



1 3 

7^ 



27r 



(50) 

(51) 
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We can check that the sound velocity obtained from the excitation spectrum 
Eq.(48) defined by 

0€(p) 

Vs = linip^o— — (52) 

gives us Vg = 2py/^ , that is exactly one order lower of the one obtained 
through the compressibility Eq.(51). 



5 Numerical Results 

We did all our calculations for a fixed p = 1. For the exact solution we 
follow the steps of Ref.|TB| solving numerically the Eqs.(5),(6) and (7). For 



the Gaussian approximation we solve the system of Eqs.(39) and (40) taking 
T = 0, and we calculate the ground state energy given by the Eq.(45). Finally 
for the pertubative results we use closed expression (49). All this results are 
shown in Fig.l, where we plot the ground state energy as a function of 7. We 
clearly see the upperbound Gaussian results and, as expected, satisfactory 
results for low 7 of both perturbative and Gaussian calculations. 

For the exact and Gaussian sound velocity we differentiate numerically 
the respective ground state curve. For the perturbative solution we have two 
possibilities one using the compressibility through Eq.(51) or in lower order 
using the excitation spectrum Eq.(48). We plot the sound velocity as a func- 
tion of 7 in Fig. 2 . We note that the Gaussian result is also an upperbound 
and for this range of 7 the results obtained from perturbative theory using 
the compressibility and the exact solution are graphically indistinguishable. 

Finally for the excitation spectrum we solve numerically (10), (11) and 
(12) for the first type of excitations and (13), (14) and (15) for the second 
type. To calculate the excitation spectrum in the Gaussian approximation 
using Eq.(39), once we have solved the system (40). For the perturbative 
theory we use Eq.(48). In Fig.3 and Fig. 4 we plot the excitation energy as 
a function of the momentum for 7 = 0.787094 and 7 = 3.07725 respectively. 
Neither perturbative nor Gaussian approximation describe the exact second 
type of excitations and the Gaussian approximation introduces a gap in the 
excitation energy that gets bigger if we go to the non-perturbative limit. 
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6 Conclusion 

In summary, we have used a one- dimensional problem to understand , in a 
clear way, the range of validity of the perturbative and the Gaussian varia- 
tional approximation in the context of interacting Bose gas. The methods 
were applied to calculated the grand-potential functional and compared to 
the exact results. In this way we avoid problems related to renormalization 



IC] and have the exact solution for the ground state energy and excitation 



spectrum [^, ^ 

We were able to see in Fig. 1 that the ground state energy for low values of 
7 both perturbative and Gaussian (non-perturbative) results are quite good. 
As 7 increases we have a region where the exact solution is in between the 
two approximations. For higher values of 7 the perturbatives results breaks 
down completelly and the Gaussian approximation still gives us an good up- 
perbound result. This shows us that for high values of 7 a self-consistent 
non-perturbative calculation should be used. For the sound velocity. Fig. 2, 
we see that both perturbative and Gaussian results are upperbound approx- 
imation to the exact result. From Figs. 3 and 4, we see a very interesting 
feature of the Gaussian approximation, namely, it produces a coherent macro- 
scopic occupation for k = 0. This occupation is artificial because as we can 
see from the exact solution there is no macroscopic occupation for A; = 0. 

Therefore, the results that are related directly from the thermodynamic 
potential, such as the ground state energy and the sound velocity, obtained 
through the compressibility, are good upperbounds. Yet the gap in the ex- 
citation spectrum, that came from the occupations, are very different from 
the exact solution. So, the fact that using the optimal determination of the 
grand potential does not assure that we will have the best variational result 



for the occupations [|T9|. Thus, for better results in other observables we 
should involve ourselves in a more difficult task of finding the best functional 
related to this specific observable with the same trial density operator of the 
form given in Eq.(20). Finally we want to point out that the formalism dev- 
eloped here alows us to go in straightfoward way to T 7^ still keeping the 
variational caracteristic of having an upperbound result |jTO| . 
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7 Figure captions 

Fig.l. The ground-state energy as a function of 7 = c/ p. The full curve 1 
gives the exact solution obtained numerically. Curve 2 show the Gaussian ap- 
proximation also obtained numerically. Curve 3 is the result of perturbation 
theory. 

Fig. 2. The sound velocity Vs-, as a function of 7 = c/ p. Curve 1 is the result 
obtained from perturbative theory using the excitation spectrum. Curve 2 
derived from the macroscopic compressibility using the Gaussian approxima- 
tion. Curves 3 and 4, that is graphically indistinguishable in this region, are 
the exact result and the perturbation solution obtained trought the macro- 
scopic compressibility. 

Fig. 3 The excitation energy as a function of the momentum for 7 = 0.787094. 
The full curves 1 and 2 corresponds to the two types of excitations (exact 
result). Curve 3 corresponds to the Gaussian spectrum and curve 4 is the 
pertubative spectrum. 

Fig. 4 The excitation energy as a function of the momentum for 7 = 3.07725. 
The full curves 1 and 2 corresponds to the two types of excitations (exact 
result). Curve 3 corresponds to the Gaussian spectrum and curve 4 is the 
pertubative spectrum. 
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